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ABSTRACT 

A series  of  measurements  of  radar  feacfcscatter  froa  rair.  have  been  cade 
at  frequencies  of  9.375,  35,  70,  and  95  GHz.  The  geometry  of  the  experi- 
ment and  the  equipment  were  chosen  such  as  to  provide  data  useful  to  the 
equipment  designer  in  the  choice  of  operating  frequency  for  his  particular 
mission.  Amplitude  statistics  for  both  linear  and  circular  polarization 
were  obtained  for  rain  rates  between  1 an/hr  and  90  en/hr.  Non-coherent 
spectral  measurements  and  correlation  properties  were  investigated  in  de- 
tail as  functions  of  frequency  and  rain  rate.  Limited  comparisons  are  made 
between  theory  and  experimental  results  and  certain  properties  of  the 
results  are  discussed  in  relation  to  the  phenomenology  of  rain  return. 


TABLE  OF  CONTENTS 


I.  INTRODUCTION 1 

A,  Background  .......  ........  1 

B.  Description  of  Radar  Field  Tests  ...........  2 

1.  Radars 2 

2.  Data  Recording  Equipment  .....  3 

3.  Calibration  Procedure.  .........  6 

4.  Radar  Reference  Targets.  ....  6 

5.  Meteorological  Instrumentation  .....  9 

6.  Measurement  Procedure.  ..............  21 

II.  DATA  ANALYSIS 13 

A.  Data  Analysis  Techniques  ........  . 13 

1.  Data  Reduction  Facility.  ...  .....  13 

2.  Data  Analysis  Procedure.  ........  13 

a.  Pulse  Height  Amplitude  Distributions  15 

b.  Frequency  Spectra.  • • , . 15 


c.  Auto-  and  Cross-Co rrelati % Functions. 


d.  Rain  Drop-Size  List  rib  uti-v.ts 


B.  Summary  of  Results  ......  .......  25 


1.  Interpretation  of  Data 


2.  Attenuation 27 

3.  Average  Return  Amplitude  .............  27 

4.  Amplitude  Fluctuations  .........  31 

5.  Power  Density  Spectrum 41 


6.  Correlation  Functions, 


III.  CONCLUSIONS  AND  RECOMMENDATIONS 59 


ijw 


LIST  OF  FIGURES 


Page 


1.  Exterior  view  of  instrument  van  showing  the  antennas 
for  the  four  radars.  ..c.  ...........  . 


5 


2.  Interior  of  instrument  van  showing  the  data  recording 
equipment.  ...................... 

3.  Block  diagram  of  experimental  equipment  set-up  for 

rain  backscatter  experiment 

4.  Field  of  view  from  radar  van  showing  radar  reference 
targets  and  rain  instrumentation  ........... 

5.  Close-up  of  4 inch  trihedral  radar  reference  target 
illustrating  deployment  of  rain  tipping  buckets.  . . . 

6.  Close-up  view  of  a rain  tipping  bucket  with  windshield 

7.  Rain  drop-size  distribution  spectrometer  deployed  '.ear 

the  4 Inch  trihedral  radar  target.  .......  .. 

8.  View  of  the  Sensor  Systems  Division’s  P2P-8/F  base.' 

data  reduction  facility.  ....  

9.  Tims  history  of  the  reco^fd  backscatter  from  rain; 
9.3?5  GHz,  33  mn/hr  rai£ - rate , W polarization  . . • . 

10.  Time  history  of  the  ^corded  backscatter  froa  rain; 

95  GHz,  33  aa/hr  rain  rate,  VY  polarization.  . . . . . 

11.  Block  diagram  of  the  equipment  setup  used  to  process 

the  rain  drsp-slze  tapes 

12.  Leasb^an  square  fit  to  radar  backscatter  data, 

V^^iarization.  

13.  Least  mean  square  fit  to  radar  backscatter  data, 

RC  polarization.  ....  ......... 

14.  Probability  density  and  cumulative  distributions 

for  backscatter  from  rain;  9.375  GHz,  5 ma/hr  rain 
rate,  W polarization 

15.  Probability  density  and  cumulative  distributions 

for  backscatter  from  rain;  35  GHz,  5 na/hr  rain 
rate,  W polarization. 

16.  Probability  density  and  cumulative  distributions 

for  backscatter  from  rain;  70  GHz,  5 mn/hr  *-ain 
rate,  W polarization.  .....  ..... 


5 

7 

8 

8 

10 

10 

14 

17 

17 

22 

29 

30 


33 


34 


35 


iii 


Page 

17.  Probability  density  and  emulative  distributions  for 
backscatter  froa  rain;  95  GHz,  5 res/ hr  rain  rate, 

VV  polarization 36 

18.  Rain  backscatter  amplitude  standard  deviatica  versus 

rain  rate  for  9.373  GHz  frequency,  W polarization  .......  37 

19.  Rain  backscatter  anplitude  standard  deviation  versus 

rain  rate  for  35  GHz  frequency,  W polarization.  ........  37 

20.  Rain  backscatter  anplitude  standard  deviation  versus 

rain  rate  for  70  GHz  frequency,  W polarization.  ........  38 

21.  Rain  backscatter  anplitude  standard  deviation  versus 

rain  rate  for  95  GHz  frequency,  W polarization 38 

22.  Rain  backscatter  anplitude  standard  deviation  versus 

rain  rate  for  9.375  Glz  frequency,  RC  polarization  .......  39 

23.  Rain  backscatter  anplitude  standard  deviation  versus 

rain  rate  for  35  GHz  frequency,  RC  polarization.  ........  39 

24.  Rain  backscatter  svlitude  standard  deviation  versus 

rain  rate  for  7 G w;  frequency,  RC  polarization.  ........  40 

25.  Rain  backscatter  a s -lit.  de  standard  deviation  versus 

rain  rate  for  95  frequency,  RC  polarization.  ........  40 

26.  Normalized  frequency  spectrun  of  rain  return  at  9.375  GHz, 

W polarization.  ........................  42 

27.  Normalized  frequency  spectrun  of  rain  return  at  35  GHz, 

W polarization.  .......................  *3 

28.  Normalized  frequency  spectrun  of  rain  return  at  70  GHz, 

W polarization.  ........................  44 

29.  Normalized  frequency  spectrun  of  rain  return  at  95  GHz, 

W polarization.  ........................  45 

30.  Normalized  frequency  spectrun  of  rain  return  at  9.375 

GHz,  RC  polarization  .............  ..  46 

31.  Normalized  frequency  spectrun  of  rain  return  at  95  GHz , 

RC  polarization.  47 

32.  Decorrelation  tine  versus  rain  rate  and  frequency  for 

rain  backscatter  49 

33.  Normalized  3uto-correlation  function  for  rain  backscatter; 

9.375  GHz  frequency,  9 se/hr  rain  rate,  W polarization.  ....  50 

iv 


Page 


Hornalized  auto-co l relation  function  for  rain 
back-scatter;  95  GHz  frequency,  9 csa/hr  rain  rate, 
W polarization.  ......  ...... 


Rain  drop-size  distributions  for  four  consecutive 
30-second  intervals  with  a light-to-aoderate  rain  rate 


Two  rain  drop-size  distributions  which  resulted  in 
different  radar  backscatter  characteristics  for  a 
heavy  rain  rate.  ................. 


Four  rain  drop-size  distributions  which  resulted  in 
different  radar  backscatter  characteristics  for  a 


exonerate  rain  rate. 


Four  rain  drop-size  distributions  which  resulted  in 
different  radar  backscatter  characteristics  for  a 
light  rain  rare 


Average  radar  backscatter  cross-section  per  unit 
volume  versu;  rain  rate;  9.375  GHz,  W polarization. 


Average  radar  backscatter  cross-section  per  unit 
voiuse  versus  rain  rate;  35  GHz,  W polarization. 

Average  radar  backscatter  cross-section  per  unit 
volune  versus  rain  rate;  70  GHz,  W polarization 

Average  radar  backscatter  cross-section  per  unit 
volune  versus  rain  rate;  95  GHz,  W polarization 


Average  radar  backscatter  crcss-section  per  unit 
volune  versus  rain  rate;  9.375  GHz,  RC  polarization. 


Average  radar  backscatter  cross-secticn  per  unit 

volune  versus  rain  rate;  35  GHz,  RC  polarization  ........  73 

Average  radar  backscatter  cross-section  per  unit 

volune  versus  rain  rate;  70  C Hz,  RC  polarization  ........  74 

Average  radar  backscatter  cross-secticn  per  unit 

volune  versus  rain  rate;  95  GHz,  RC  polarization  ........  75 

Amplitude  of  the  frequency  spectrum  of  a logarithmic 
reciever  for  rain  backscatter;  5 na/hr  rain  rate, 

9.375  GHz  frequency,  and  W polarization  78 


Aoplitude  of  the  frequency  spectrua  of  a logarithmic 
receiver  for  rain  backscatter;  5 an/hr  rain  rate, 

35  GHz  frequency,  and  W polarization...  ...... 


v 


1 


B3. 

B4. 

B5. 

B6. 

37. 

38. 
B9. 
BIO. 
Bll. 
B12. 
B13. 
B14. 
E15. 


paSc'. 


Amplitude  of  Che  frequency  speccrun  of  a logarithnic 
receiver  for  rain  backscatter;  5 crs/hr  rain  race, 

70  GHz  frequency,  and  VY  polarization.  79 

Amplitude  of  Che  frequency  speccrun  of  a logarithmic 
receiver  for  rain  backscatzer;  5 m/hr  rain  race, 

95  GHz  frequency,  and  W polarization 79 

Anplitude  of  the  frequency  spectrun  of  a logarithnic 
receiver  for  rain  barkscatter;  23  cars/hr  rain  rate, 

9.375  GHz  frequency,  and  W polarization  ....,  8° 

Amplitude  of  the  frequency  spectrua  of  a logarithmic 
receiver  for  rain  backscatter;  23  m/hr  rain  rate, 

35  GEiz  frequency,  and  W polarization.  .............  80 

Anplitude  of  the  frequency  spectrua  of  a logarithmic 
receiver  fer  rain  backscatter;  23  nrs/hr  rain  rate, 

70  Giz  frequency,  and  W polarization.  

Anplitude  of  the  frequency  spectrun  of  a logarithmic 
receiver  for  rain  backscatter;  23  na/hr  rain  rate, 

95  Giz  frequency,  and  W polarization.  ........  81 

Anplitude  of  the  frequency  spectrun  of  a logarithmic 
receiver  for  rain  backscatter;  38  sc/hr  rair*  rate, 

9.375  GHz  frequency,  and  W polarization 32 

Amplitude  of  the  frequency  spectrun  ef  a logarithmic 
receiver  for  rain  backscatter;  38  m/hr  rain  rate, 

35  Giz  frequency,  and  VV  polarization. 82 

Amplitude  of  the  frequency  spectrun  of  a logarithmic 
receiver  for  rain  backscatter;  38  oes/hr  rain  rate, 

70  Giz  frequency,  and  W polarization.  .............  83 

Amplitude  of  the  frequency  spectrua  of  a logarithmic 
receiver  for  rain  backscatter;  38  szs/hr  rain  rate, 

95  Giz  frequency,  and  W polarization.  ........  83 

Anplitude  of  the  frequency  spectrua  of  a logarithmic 
receiver  for  rain  backscatter;  70  an/hr  rain  rate, 

9.375  GHz  frequency,  and  W polarization.  8* 

Anplitude  of  the  frequency  spectrua  of  a logarithnic 
receiver  fer  rain  backscatter;  70  sw/hr  rain  rate, 

35  GHz  frequency,  and  W polarization . ...  84 

Anplitude  of  the  frequency  speccrun  of  a logarithnic 
receiver  for  rain  backscatter;  70  Ro/hr  rain  rate, 

70  Giz  frequency,  and  W polarization .85 


vi 


516.  Amplitude  of  the  frequency  spectrins  of  a logarithnic 
receiver  for  rain  backscatter;  70  na/hr  rain  rate, 

95  GHz  frequency,  znd  W polarization 

317.  Amplitude  of  the  frequency  s pec t run  of  a logarithnic 
receiver  for  rain  backscatter;  100  m/hr  rain  rate, 

9.375  GHz  frequency,  and  W polarization  . ...» 

BIS.  Amplitude  of  the  frequency  spectres  of  a logarithnic 
receiver  for  rain  backscatter;  100/zc/hr  rain  rate. 
35  GHz  frequency,  and  W polarization.  ....... 

B19.  Amplitude  of  the  frequency  spectrins  of  a logarithnic 
receiver  for  rain  backscatter;  100  rs/hr  rain  rate, 
70  GHz  frequency,  and  W polarization.  ....... 

320.  Amplitude  of  the  frequency  spectrus  of  a logarithnic 
receiver  for  rain  backscatter;  100  m/hr  rain  rate, 
95  GHz  frequency,  and  W polarization.  ....... 

321.  Amplitude  of  the  frequency  spectrins  of  a logarithnic 
receiver  for  rain  backscatter;  3.1  m/hr  rain  rate 

9.375  GHz  frequency,  and  PjC  polarization  ...... 

522..  Amplitude  of  the  frequency  spectrus  of  a logarithnic 
receiver  for  rain  backscatter;  3.1  ess/fcr  rain  rate, 
35  C !z  frequency,  and  RC  polarization 

B23.  Amplitude  of  the  frequency  spectrum  of  a logarithnic 
receiver  for  rain  backscatter;  3.1  an /hr  rain  rate, 
70  GHz  frequency,  and  ?-C  polarization . ....... 

B24.  Amplitude  of  the  frequency  spectrins  of  a logarithnic 
receiver  for  rain  bacsscatter;  3.1  sn/hr  rain  rate, 
95  GHz  frequency,  and  RC  polarization.  ....... 

322  Amplitude  of  the  frequency  spectrum  of  a logarithnic 
-eceiver  for  rain  backscatter;  32  m/hr  rain  rate. 

9.375  GHz  frequency,  and  RC  polarization  ...... 

326.  Amplitude  of  the  frequency  spectrus  of  a logarithmic 
receiver  for  rain  backscatter;  32  m/hr  rain  rare. 

35  GHz  frequency,  and  RC  polarization 

327.  Amplitude  or  the  frequency  spectrum  cf  a logarithnic 
receiver  for  rain  backscatter:  32  m/hr  rain  rate, 

70  GHz  frequency,  and  RC  polarization 

B28.  Asplitude  of  the  frequency  spectrus  cf  a logarithnic 
receiver  for  rain  backscatter;  32  sxs/lir  rain  rate, 

95  GHz  frequency,  and  RC  polarization.  ....... 


Page 

85 

86 
86 

87 
S' 

88 
88 
89 
89 
SO 
"•0 
91 
91 


vii 


Page 


B29.  Amplitude  of  the  frequency  spectrum  of  a logarithmic 
receiver  for  rain  backscatter;  91  mm/hr  rain  rate, 

9.375  GHz  frequency,  and  RC  polarization. 92 

B30.  Amplitude  of  the  frequency  spectrum  of  a logarithmic 
receiver  for  rain  backscatter;  91  mm/hr  rain  rate, 

35  GHz  frequency,  and  RC  polarization 92 

B31.  Amplitude  of  the  frequency  spectrum  of  a logarithmic 
receiver  for  rain  backscatter;  91  mm/hr  rain  rate, 

70  GHz  frequency,  and  RC  polarization  .....  93 

B32.  Amplitude  of  the  frequency  spectrum  of  a logarithmic 
receiver  for  rain  backscatter;  91  mm/hr  rain  rate, 

95  GHz  frequency,  and  RC  polarization  .....  ....  93 


LIST  OF  TABLES 


Page 

Table  1.  Radar  Parameters  for  Rain  Experiment  [3],  . . . 4 

Table  2.  Attenuation  Due  to  Rain  ............ 28 


ix 

g - ■ - - 


I.  INTRODUCTION 


This  report  summarizes  the  results  of  an  analysis  of  measurements  on 
the  radar  backscatter  from  rain  at  millimeter  frequencies  which  was  con- 
ducted at  McCoy  AFB,  Florida,  in  August  and  September  1973  by  the  U.  S. 

Army  Ballistic  Research  Laboratories  (BRL) , Aberdeen,  Md.  Emphasis  in  the 
report  is  placed  on  those  characteristics  of  the  radar  backscatter  which 
could  be  determined  from  data  recorded  on  magnetic  tape,  as  opposed  to  a 
detailed  description  of  the  actual  radar  tests  or  analysis  of  photographic 
and  video  tape  data  which  were  also  gathered.  These  parameters  will  be 
discussed  elsewhere  in  a separate  report  to  be  published  by  BRL. 

A.  Background 

Under  Contract  DAA25-73-C-0256  with  the  U.  S.  Army,  Frankford  Arsenal, 
the  Engineering  Experiment  Station  (EES)  at  Georgia  Tech  undertook  to 
develop  suitable  mathematical  models  to  allow  computer  simulation  of  fire- 
control  radar  systems  in  the  millimeter  region.  As  a first  step  in  this 
program,  a literature  search  was  conducted  to  determine  what  data  were 
available  at  millimeter  frequencies  on  clutter  and  target  characteristics 
which  would  affect  system  performance.  A summary  of  these  data  and  the 
empirical  models  which  were  developed  were  presented  in  Technical  Report  1 
on  the  contract  (1], 

The  results  of  that  study  brcught  to  light  a number  of  deficiencies 
in  the  available  data.  It  was  determined  that:  (1)  little  or  no  data  are 
available  on  the  characteristics  of  ground  clutter  above  10  GHz,  particularly 
below  1Q°  grazing  angle;  (2)  almost  no  data  were  available  on  the  effects 
of  atmospheric  conditions  such  as  precipitation  in  the  millimeter  region; 
and  (3)  there  was  a scarcity  of  data  on  the  radar  cross-section  (RCS) 
properties  of  hard  targets  at  these  frequencies.  EES  then  proposed  to 
Frankford  Arsenal  a program  of  investigation  to  fill  in  these  gaps  in  the 
available  data  to  be  followed  by  an  analysis  program  aimed  at  tying  the 
data  results  together  in  a unified  model  for  millimeter  radar  systems. 

Frankford  Arsenal,  in  response  to  this  proposal,  undertook  to  fund 
studies  of  the  backscatter  from  precipitation  and  from  land  clutter  at 
millimeter  frequencies. 


Ballistic  Research  Laboratories  was  funded  by  Frankford  to  conduct 
an  experiment  on  backscatter  from  rain  at  10  GHz,  35  GHz,  70  GHz,  and 
95  GHz,  and  EES  was  tasked,  in  a modification  to  the  original  Army  con- 
tract, to  assist  in  the  experiment  planning  and  to  act  as  on-site  observers 
during  the  actual  tests.  In  addition,  under  this  program,  EES  constructed 
and  delivered  to  BRL  a range-gated  boxcar  sampler  unit  which  would  allow 
recording  of  backscatter  d-vC  on  magnetic  tape  during  the  tests.  The 
services  and  activities  performed  were  detailed  in  a letter  report  at 
the  end  of  the  field  tests  12 J. 

After  the  termination  of  the  radar  fie Ld  tests  on  rain  backscatter, 
Frankford  Arsenal  provided  additional  funds  to  EES  for  the  purpose  of 
reducing  and  analyzing  the  data  contained  on  the  magnetic  tapes  recorded 
during  the  tests.  The  goals  of  this  analysis  program  were  to  (1)  determine 
average  backscatter  parameters  in  terms  of  rain  rate  and  drop-size  distri- 
butions as  a supplement  and  cross  check  to  the  results,  obtained  by  BRL 
and  (2)  determine  spectral  auto-correlaticn  and  o.her  properties  of  the 
backscatter  which  were  not  obtained  during  the  i.  Ltial  analysis  which 
used  primarily  photographic  methods.  This  report  summarizes  the  results 
of  the  magnetic  tape  analysis  program. 

B.  Description  of  Radar  Field  Tests  1 31 

The  radar  field  tests  were  conducted  at  McCoy  Air  Force  3ase,  Orlando, 
Florida,  during  August  and  Xeptcraber  1973.  This  ->ii-?  was  selected  because 
of  the  afternoon  rains  that  occur  almost  daily  in  this  area  during  the 
summer  months.  During  the  tests,  it  rained  25  tines,  the  rainfall  rates 
varying  from  light  drizzles  up  to  very  heavy  rainfalls  of  150  ure/hr. 

For  this  experiment,  an  instrument  van  containing  the  four  radars 
and  the  data  recording  equipment  was  positioned  at  one  end  of  a.  large 
field,  and  the  radar  targets  and  rain  instrumentation  were  placed  near 
the  other  end.  Data  were  taken  whenever  a rain  celi  passed  through  the 
region  around  the  targets. 

1.  Radars 

Four  radars  .'ere  used  during  the  rain  experiments  with  operating 
frequencies  of  9.375  GHz,  35  GHz,  70  GHz , and  95  GHz,  The  basic  radar 


parameters  are  given  in  Table  1 [3],  Logarithmic  receivers  were  provided 
for  all  of  the  radars,  and  0 to  50-dB  variable  waveguide  attenuators  were 
installed  in  the  rf  paths  to  the  receivers  for  calibration  purposes. 

Parabolic  reflector  antennas  with  Cassegrain  feeds  provided  pencil 
beams  for  all  the  radars.  Vertical,  right  circular,  and  left  circular 
polarizations  were  available.  The  antennas  were  attached  to  mounts,  which 
could  be  manually  positioned  in  azimuth  and  elevation,  and  boresighted 
telescopes  were  provided  for  each  mount.  The  radars  and  antennas  were 
housed  in  an  air-conditioned  instrument  van  with  a 1/8-inch  thick 
Polycast  window  along  one  side  to  serve  as  a radome.  Figure  1 gives  an 
exterior  view  of  the  instrument  van  showing  the  antennas.  The  antennas, 
as  shown  left  to  right,  are  for  the  S.375  GHz,  35  GHz,  70  GHz,  and  95  <3iz 
radars  respectively. 

2.  Data  Recording  Equipment 

During  the  tests,  data  were  recorded  on  the  rain  backscatter  in 
several  ways  in  order  to  obtain  as  much  information  as  possible.  The 
received  signals  from  each  radar  were  displayed  on  an  A-scope  and 
Polaroid  photographs  were  made  at  regular  intervals.  In  addition,  video 
recordings  were  made  on  the  A-scope  displays,  changing  the  video  displayed 
at  10-second  intervals.  Finally  the  videos  from  the  four  radars  were 
sampled  by  a four-channel,  range-gated  boxcar  unit  which  stretched  the 
return  from  a narrow  gate  so  as  to  allow  recording  on  lew-bandwidth 
devices  such  as  a strip  chart  recorder  and  an  FM  magnetic  tape  recorder. 

In  order  to  allow  simultaneous  data  taking,  all  of  the  radar  modu- 
lators were  synchronized  to  a common  trigger  source.  To  accomplish  this, 
the  1600-Hz  prf  from  the  95-GHz  radar  timer  was  used  to  trigger  the 
9.375-GHz  and  70-GHz  radars  directly,  and  was  doubled  to  3200  Hz  to  trigger 
the  35-GHz  radar.  Also,  the  four  range  gates  i rf  the  boxcar  sampler  were 
slaved  together  and  controlled  by  a single  range  delay  adjustment,  so  that 
simultaneous  data  could  be  recorded  at  a given  range  at  all  four  fre- 
quencies. 

Figure  2 gives  a view  of  the  inside  of  the  instrumentation  van  showing 
the  A-scopes,  video  monitor,  boxcar  sandier,  FM  magnetic  tape  recorder. 
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time-code  generator,  and  strip  chart  recorder,  while  Figure  3 shows 
a block  diagram  of  the  data  instrumentation  setup. 

The  strip  chart  recorder  served  as  a monitor  on  the  radar  performance, 
and  was  also  used  during  the  tests  to  record  the  bucket  tips  from  the  rain 
tipping  buckets  for  later  analysis.  During  the  data  reduction  at  a 
later  time,  the  charts  were  particularly  helpful  in  identifying  data  on 
the  FM  magnetic  tapes,  and  in  providing  supplementary  information  in 
addition  to  the  data  sheets. 

3.  Calibration  Procedure 

The  variable  attenuators  installed  in  the  t-f  signal  path  of  the  radars 
were  used  to  calibrate  all  the  data.  The  procedure  consisted  of  bore- 
sighting all  four  radars  on  to  one  of  the  calibrated  radar  targets  (usually 
a 4-inch  tricomer  reflector)  and  recording  the  return  from  the  target 
as  the  attenuators  were  varied  from  0 to  50  dB  in  5 dB  steps.  In  the  case 
of  the  35  GHz  radar,  the  receiver  was  saturated  by  the  return  from  the 
comer,  so  that  it  is  necessary  to  insert  an  extra  30  dB  pad  in  series 
with  the  variable  attenuator  and  repeat  the  above  procedure.  The  cali- 
brations were,  of  course,  performed  before  and/or  after  a rain  store  when 
no  rain  was  falling  in  order  to  eliminate  the  effects  of  rain  attenuation 
on  the  calibration.  All  the  attenuators  were  calibrated  against  laboratory 
standards  after  the  experiment. 

4.  Radar  Reference  Targets 

A 4-inch  trihedral  comer  reflector  was  mounted  on  a 32-foot-high 
fiberglass  mast  in  the  field  in  front  of  the  radars  at  a range  of  1478  feet 
for  use  as  a primary  reference  target.  Figure  4 is  a view  of  the  field 
from  the  radar  van,  showing  the  reference  targets  and  the  rain  instrumen- 
tation. The  primary  reference  target  was  chosen  to  be  large  enough  to 
provide  sufficient  return  to  overlde  ground  clutter  and  backscatter  from 
heavy  rain,  but  not  so  large  that  the  signal  could  not  be  attenuated  down 
to  the  receiver  noise  level  with  available  attenuators.  The  aluminum 
trihedral  was  accurately  machined,  and  had  a radar  cross-section  of 
0.43,  7,  16,  and  45  square  meters  at  9.375,  35,  70,  and  95  Qlz,  respectively. 
The  trihedral  was  calibrated  against  a 10-inch-diaaeter  aluminum  sphere 
suspended  from  a rope  strung  between  two  poles.  Figure  5 shows  a close-up 
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of  che  trihedral  in  the  field  viewed  froa  the  radar  van.  1-93  h2ig‘nt  of  the 
reference  trihedral  was  a compromise  intended  to  fce  groat  enough  to 
minimize  ground  clutter  and  multipath  effects  but  not  so  ,reat  as  to  be 
unstable  in  orientation  during  a heavy  rain  and  wind  storm.  The  range  of 
the  reference  target  from  the  radars  was  determined  by  three  criteria!  (1) 
having  the  range  as  short  as  possible  to  maximize  the  elevation  angle  to 
the  trihedral,  thereby  minimizing  ground  clutter,  multipath,  and  vegetation 
interference  (the  latter  was  initially  a problem  with  6-foot  high  weeds 
in  the  field  which  had  to  be  cut  by  hand) ; (2)  having  the  range  as  great 
as  possible  to  simulate  long-range  radar  operation  and  still  have  measurable 
rain  backscatter  during  all  rainfall  intensities;  and  (3)  having  the 
target  beyond  the  near-field  of  the  antennas,  (h'ote:  the  near-field  criterion 
was  not  achieved  for  the  95-GHz  radar.)  The  target  height  of  32  feet  and 
range  of  1473  feet  resulted  in  no  discernible  ground  clutter  at  35,  70,  and 
95  GHz  at  the  target  range;  the  ground  clutter  at  9.375  Ciz  was  45  dB  down 
from  the  target  return.  Other  refernce  targets  which  were  also  used  from 
time  to  time  included  a 6-inch  tricorner  and  a 2-inch  dihedral. 

5.  Meteorological  Instrumentation 

Jleceorologicai  conditions  which  were  recorded  during  the  tests  include 
rainfall  rate  (from  the  tipping  buckets),  rainfall  rate  and  drop-size  dis- 
tributions (from  the  rain  spectrometers) , temperature,  atmospheric  pressure, 
relative  humidity,  and  wind  speed  and  direction.  The  rain  tipping  buckets 
and  spectrometers  were  located  near  the  radar  targets,  while  all  other 
weather  instrumentation  was  located  at  the  radar  van. 

The  basic  measure  of  rainfall  rate  ves  obtained  from  the  rain  tipping 
buckets.  Fig-ire  6 gives  a close-up  view  of  a tipping  bucket  with  wind- 
shield. (TVo  of  the  buckets  had  windshields  and  one  did  not.)  A wind- 
shield consisted  of  a ring  of  hinged  vanes  intended  to  breakup  cross 
winds.  Each  tipping  bucket  had  a 7-inch  diameter  orifice  and  provided 
a switch  closure  for  each  .01  inch  (0.254  ejk)  of  rainfall.  As  previously 
mentioned,  the  switch  closures  from  each  bucket  were  recorded  as  pulses 
on  the  paper  chart  recorder.  It  has  since  been  discovered  that  the 
accuracy  of  the  tipping  buckets  in  determining  rain  rate  is  in  doubt 
below  10  aa/hr  rain  rates  and  above  60  na/hr  rain  rates. 
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The  rain  drop-size  spectrometers  were  located  near  the  tipping  buckets  , 
and  were  intended  to  provide  data  on  the  rain  particle  sizes  and  to  serve 
as  a check  on  the  rain  races  indicated  by  the  tipping  buckets.  Figure  7 
gives  a close-up  view  of  one  of  the  spectrometers.  These  units  were 
furnished  by  the  Illinois  State  Water  Survey,  Urbana,  Illinois,  and  con- 
sisted of  a piezoelectric  sensor  mounted  under  a plastic  cylinder,  2-3/4  inches 
long.  The  sensor  generated  a voltage  pulse  for  each  raindrop  that  impacted 
on  the  top  surface  of  the  plastic  cylinder,  with  an  amplitude  that  was  pro- 
portional to  the  drop  size.  The  sensor  also  contained  a preamplifier  and 
an  electronic  switch  so  that  a calibration  signal  could  be  introduced 
periodically. 

A cassette  tape  recorder  activated  by  the  first  raindrop  to  hit  the 
spectrometer,  was  used  to  record  the  voltage  pulses  along  with  a time-code 
signal.  A field  calibration  of  the  spectrometer  was  achieved  by  dropping 
a plastic  bead  from  . known  height  onto  the  spectroneter  and  recording  the 
resulting  pulse.  Correlation  between  the  bead-drop  calibration  and  water 
drop-size  was  done  in  a water  drop  tower  test  by  the  Illinois  Water  Survey 
Laboratory,  who  also  reduced  a portion  of  the  drop-size  tapes  following  the 
field  operations.  An  important  limitation  of  Che  spectrometers  was  their 
inability  to  measure  drops  below  0.5  na  diameter. 

The  other  meteorological  instruments  were  located  at  the  instrument 
van,  and  readings  such  as  temperature,  barometric  pressure,  and  wind  speed 
and  direction  were  cade  at  intervals  during  a rain  storn. 

6.  Measurement  Procedure 

The  measurement  procedure  followed  in  taking  rain  backs catter  data 
typically  consisted  of  the  following  steps: 

(1)  At  the  onset  of  rain  (which  was  usually  detected  by  the  70  or 
95  GHz  radars  well  iu  advance  of  the  tipping  cf  the  buckets),  all  the 
recording  equipment  was  turned  on. 

(2)  A-scope  photographs  were  taken  at  intervals  during  the  rain  storn, 
as  considered  necessary  to  record  the  complete  range  of  rain  intensities 
being  monitored  via  the  bucket  tips  recorded  on  the  paper  chart  recorder. 
Typically,  A-scope  photographs  were  £3ken  at  intervals  of  between  one  and 
five  ainutes.  Photos  were  s>ade  for  all  three  available  polarizations 
during  each  interval. 


11 


(3)  Ac  5-  co  xO-cdRute  intervals  the  range  gate  controlling  the 
pulse  sampler  boxcai  was  changed  to  different  settings,  thus  providing 
a magnetic  Cape  and  paper  chart  record  of  the  amplitude  of  the  back- 
scatter  versus  rain  rate  for  different  ranges » 

(4)  The  output  oi  each  radar  was  switched  sequentially  at  about 
10-seccnd  intervals  to  the  A-scepe  being  recorded  on  video  tape. 

(5)  Temperature,  pressure,  relative  humidity,  and  wind  speed  and 
direction  were  measured  at  intervals  during  the  rain  store, 

(6)  Calibration  of  all  recorders  was  performed  before  and/or  after 
the  rain,  the  radars  being  boresighted  on  the  reference  target  and  the 
rf  attenuation  varied  in  5-dB  steps. 
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11.  DATA  ANALYSIS 
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A.  Data  Analysis  Techniques 

The  data  analysis  performed  under  this  contract  consisted  in  reducing 
the  data  recorded  on  magnetic  tapes  during  the  field  operations  at  McCoy 
AFB.  Data  from  .other  sources,  such  as  Polaroid  photographs  of  A-scope 
displays,  video  Recordings,  and  strip  chart  plavouts,  served  to  provide 
background  information  but  were  not  analyzed  directly. 

1,  Data-P.cduction  Facility 

The  FDF-8/F  based  data- reduct ion  facility  of  the  Sensor  Systens 
Division  of  EKS  was  used  to  process  all  the  tapes  which  resulted  from 
the  McCoy  AFB  field  tests.  Figure  8 gives  a view  of  the  basic  computer 
components.  These  include:  (1)  An  analog  signal-conditioner  unit  which 
provides  variable  gain  an  i offset  to  allow  the  interface  of  varied  types 
of  signals  to  the  data-reduction  facility,  (2)  A Fabritek  Series  1070 
instrument  computer  which  serves  as  an  A/D  and  D/A  interface,  and  also 
computes  real-time  pulse-height  distributions  and  cress-correlation 
functions.  The  D/A  output  from  the  Fabritek  computer  can  be  dispi  yed  on 
a CRT  display  or  be  plotted  on  an  x-y  plotter.  (3)  A PDP-8/F  computer 
which  can  exchange  information  directly  with  the  Fabritek  computer.  (4) 

A teletype,  and  (5)  a Sykes  Compucorp  Digital  Cassette  Recorder  for  pro- 
gram development  and  storage. 

Tne  PDP-8/F  contains  12K  of  memory,  of  which  8K  is  magnetic  core. 

An  extended  version  of  F0CALTV,  has  been  developed  for  use  with  the 
PDP-8/F  and  is  designated  FOCL/F  [4],  This  language  is  interactive  and 
greatly  facilitates  program  correction  and  modification.  Also  available 
is  a machine  language  software  package  for  calculating  fast  Fourier  trans- 
forms (FFT),  and  a set  of  software  coianands  for  Fabritek  control.  These 
two  machine  language  software  packages  along  with  the  extended  FOCAL^ 
software  make  this  system  a very  powerful  and  flexible  data-reduction 
tool. 

2.  Data  Analysis  Procedure 

The  types  of  results  which  can  be  obtained  from  the  data  facility 
include:  (1)  pulse-height  amplitude  distributions,  (2)  frequency  spectra. 
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of  the  Sensor  Systems  Division's  PDP-8/F  based 
reduction  facility. 


and  (3)  auto-  and  cross-correlation  functions*  The  methods  for  obtaining 
these  three  classes  of  results  are  sufficiently  different  as  to  require 
entirely  different  FOCAL_,  programs  for  their  calculation.  Taken  as  a 
whole,  these  three  results  are  considered  to  describe  the  recorded  data. 

a.  Pulse-Height  Amplitude  Distributions 

Pulse-height  amplitude  distributions  calculated  by  the  data  facility 
are  displayed  in  two  forms:  as  probability  density  plots  and  as  probability 
distribution  functions.  The  probability  density  plots  are  generated 
from  data  time  histories  using  a Fabritek  Type  SW-75  PHA  plug-in  as  a 
preprocessor.  This  plug-in  samples  the  input  analog  time  history,  A/D 
converts  the  samples,  determines  into  which  of  1024  amplitude  bins  the 
sample  belongs,  and  increments  a stored  variable  corresponding  to  the 
number  of  samples  which  have  fallen  in  that  amplitude  windcw.  When 
repeated  a large  number  of  timet,  this  process  generates  an  unnormalized 
density,  which  is  then  divided  by  the  total  number  of  samples  to  achieve 
a normalized  probability  density  function. 

The  voltage  amplitudes  are  calibrated  by  inputting  a stepped  cali- 
bration signal  to  the  PHA  program.  The  peak  of  the  distribution  for  each 
voltage  step  in  the  calibration  is  then  assigned  the  dB  value  corresponding 
to  that  calibration  step.  The  PHA  program  in  the  PDP-8/F  then  does  a 
cubic  fit  to  the  calibration  and  generates  a table  relating  dB  value 
and  amplitude  bin  number.  The  cubic  fit  program  was  developed  to 
"linearize"  nonuniform  calibration  steps  so  that  the  output  density  func- 
tions can  be  plotted  on  a linear  scale. 

The  probability  distributions  are  calculated  by  a point-by-point 
numerical  integration  of  the  probability  density  functions.  The  func- 
tional values  of  these  distributions  are  then  multiplied  by  a nonlinear 
transfer  function  so  that  they  can  be  plotted  on  probability  paper. 

The  resultant  probability  distributions  are  useful  in  determining 
the  median  values  of  the  distributions  and  also  their  shapes.  In  addition, 
for  certain  classes  of  functions,  the  average  values  can  be  determined 
from  the  distributions  using  a simple  formula  [5]. 

b.  Frequency  Spectra 

The  frequency  spectra  data-reduction  program  uses  the  fast  Fourier 
transform  subroutines  available  for  the  PDF-8/F  to  transform  input  time 
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histories  to  the  frequency  domain.  The  program  has  several  options 
which  allow  the  spectral  data  to  be  presented  in  several  different  forms 
including:  (1)  voltage  amplitude  in  dB  versus  frequency,  (2)  normalized 
voltage  amplitude  in  dB  versus  frequency,  and  (3)  normalized  power 
spectral  density  versus  frequency. 

The  voltage  amplitude  program  in  the  PDP-8/F  computer  processes  an 
input  time  history  which  has  been  sampled,  A/D  converted,  and  stored  by 
the  Fabritek.  The  program  computes  the  fast  Fourier  transform  for  that 
time  history,  and  calculates  the  square  roof  of  the  sun  of  the  squares 
of  the  real  and  imaginary  parts  of  each  element  in  the  transform.  Since 
the  time  history  represents  the  voltage  out  of  a logarithmic  receiver, 
the  amplitude  of  the  spectrum  corresponding  to  the  time  history  is 
measured  in  dB  relative  to  a milliwatt  (dBm).  Thus,  to  calibrate  the 
spectral  amplitude  in  dB,  a calibration  is  stored  in  memory  which  relates 
dB  values  to  input  voltage  amplitudes,  and  each  time  history  is  converted 
to  dB  values  as  it  is  read  in.  Since  the  Fourier  transform  is  a linear 
process,  the  spectrum  amplitude  will  be  proportional  to  dB  if  the  tine 
history  is  calibrated  in  dBm,  and  the  proportionality  constant  is  set 
equal  to  one  in  the  program. 

' Figures  9 and  10  give  typical  time  histories  for  9.375  CK2  and  95 

GHz.  The  vertical  scales  are  calibrated  in  units  of  dB  relative  to 
2 3 

1 m /m  . (Received  power  can  be  related  to  cross-section  per  unit 
volume  if  the  appropriate  radar  constants  are  known.)  A time  history 
is  limited  in  duration  to  the  number  of-  Fabritek  memory  bits  times  the 
sample  period.  For  this  case,  the  sample  rate  equals  the  prf  used  in 
taking  the  data,  so  that  a time  history  is  given  by:  1/1620  x 1024  = 
0.632  seconds.  To  achieve  the  equivalent  of  a longer  time  history,  the 
Fourier  transforms  of  8 adjacent  time  histories  are  averaged  together. 

In  calculating  the  FFT,  the  dc  term  is  set  to  zero  because  its 
amplitude  is  normally  so  much  larger  than  the  rest  of  the  spectrum  that 
dynamic  range  problems  are  encountered  in  the  computer.  Thus,  the 
zero  frequency  point  is  zeroed  in  all  the  plots.  However,  this  dc  point, 
can  still  be  determined  independently  from  the  amplitude  distributions. 

A second  method  of  displaying  the  log  amplitude  spectrum  is  to 
normalize  the  spectrum  by  dividing  all  the  elements  by  the  peak  element 
value.  The  spectrum  is  then  plotted  as  dB  relative  to  the  peak  voltage 
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amplitude,  which  is  normally  the  lowest  frequency  po.  ,.t.  This  type  of 
plot  is  very  useful  for  comparing  frequency  .-olloff  characteristics  of 
spectral  plots  with  different  amplitudes. 

Although  the  data  recorded  on  the  magnetic  tapes  represents  the 
voltage  output  from  a logarithmic  amplifier,  it  was  deemed  desirable 
to  be  able  to  plot  an  equivalent  spectrum  for  a linear  receiver.  There- 
fore a modification  was  developed  to  the  spectrum  program  which  allowed 
the  "delogging”  of  the  input  data  prior  to  calculation  of  the  Fourier 
transform.  When  normalized  this  program  results  in  the  calculation  of 
the  normalized  power  spectrum.  For  certain  classes  of  functions,  to 
which  all  of  the  data  processed  here  belong,  tbe  power  spectrum  is 
equivilent  to  the  power  spectral  density.  The  unnormalized  calibrated 
power  spectrum  cannot  be  determined  because  the  dc  term  is  thrown  away, 
as  previously  described. 

Another  program  modification  was  developed  to  de3l  with  both  60-Hz 
and  400-Hz  interference  signals  in  the  recorded  data.  The  400-Hz  (approxi- 
mately) component  was  often  so  much  larger  than  tne  backscatter  data  that  the 
maximum  allowed  dynamic  range  for  the  computer  variables  was  exceeded.  To 
deal  with  this  problem,  the  tine  history  was  passed  through  a six-pole  Butter- 
worth  notch  filter  with  a center  frequency  of  400  Hz  prior  to  being  A/J)  con- 
verted. An  algorithm  for  the  inverse  of  the  amplitude  versus  frequency- 
response  of  the  filter  was  developed  which,  when  multiplied  by  the 
spectrum  of  the  filtered  data,  resulted  in  a spectrum  equivalent  to  the 
uafiltered  spectrum.  The  notch  filter  decreased  the  size  of  the  400-Hz 
component  relative  to  the  other  spectral  components,  thus  solving  the 
dynamic  range  problem.  The  value  of  the  400-Hz  component  was  then  set  to 
zero  before  the  inverse  filter  algorithm  was  applied,  yielding  a net  result 
of  a spectrum  equivalent  to  the  unfiltered  spectrum  except  for  the  absence 
of  the  400-Hz  component.  The  400-Hz  signal  varied  in  frequency  as  seen  in 
Appendix  B. 

C.  Auto-  and  Cross-Correlation  Functions 


The  auto-correlation  functions  and  cross-correlations  functions  for 
input  time  histories  are  determined  by  two  different  methods  in  the  data 
analysis.  Auto-correlation  functions  are  computed  from  the  inverse 
transform  of  the  magnitude  of  the  Fourier  transform  squared  in  the  PDP-8/F 
computer,  while  cross-correlation  functions  are  measured  using  an  SD-75 
plug-in  in  the  Fabritek  computer.  The  SD-75  plug-in  samples  two  separate 


input  signals,  A/D  converts  them,  and  calculates  partial  products  between 
samples  as  a function  of  time  to  generate  the  correlation  function  between 
the  signals.  This  computation  is  given  by 


*AB(niT) 


iV  , 

k JLj  ' 


A(t^)  + n£t);  n * 0,  l,2,...,k-i 


where,  k * number  of  discrete  points  for  which  the  correlation 
function  is  to  be  determined, 
kAt  3 total  sample  length, 

A(t.J  * value  of  first  sampled  function  at  time  t^ , 

X X 

B(tt.)  » value  of  second  sampled  function  at  time  t^. 


This  method  could  have  been  used  to  generate  autocorrelation 
functions  by  connecting  the  same  signal  to  both  plug-in  inputs,  but  this 
was  not  done  because  of  several  limitations  in  the  plug-in.  These  include: 
(1)  There  is  co  provision  to  calibrate  the  twe*  input  signals,  making 
set-up  and  calibration  of  the  auto-correlation  function  difficult.  (2) 

Due  to  the  liaised  memory  available  in  the  Fabritek,  the  lowest  frequency 
that  can  be  accurately  measured  is  approximately  one-tenth  the  total 
sampling  period  while  the  highest  frequency  is,  from  the  sailing  theorem, 
about  one-half  the  sampling  rate.  Thus  to  measure  the  auto-correlation 
function  of  3 signal  with  800-Hz  bandwidth,  the  sample  rate  would  be 
1600  Hz  and  the  lowest  frequency  would  be: 


„ m Sample  rate  m 1600  H» 

'low  Total  Memory  1024  Memory 

bits 


1.56  Hz. 


(3)  The  finite  computation  time  limits  the  sampling  frequency  since  com- 
putations are  performed  in  real  time. 

Utilizing  the  fast  Fourier  transform  in  the  PDP-8/F  to  perform  auto- 
correlation calculations  solves  the  problems  inherent  in  using  the  SD-73 
plug-in  in  the  Fabritek,  except  for  the  lower  frequency  limitation.  Since 
the  Fabritek  memory  is  used  to  store  the  time  histories,  the  same  low- 
frequency  limitations  apply.  However,  a partial  solution  to  this  problem 
is  to  calculate  the  auto-correlation  function  using  a lower  sample  rate 


to  pick-up  lower  frequencies  at  the  expense  of  higher  frequencies.  A 
family  of  curves  can  thus  be  generated  at  different  sampling  rates  which 
describe  the  auto-correlation  function  over  any  desired  band  of  frequencies. 

In  the  case  of  cross-correlation  function  computation,  the  FFT  pro- 
gram in  the  PDP-8/F  computer  is  not  suitable,  and  the  SD-75  plug-in  in 
the  Fabritek  oust  be  used  even  with  its  limitations.  Calibration  can  be 
obtained  with  the  plug-in  by  using  signals  with  known  correlation  func- 
tions, such  as,  for  example  a sinewave  or  Gaussian  noise,  to  calibrate 
the  results  from  an  unknown  signal. 

d.  Rain  Dron-Size  Distributions 

The  rair.  drop-size  distribution  analysis  was  carried  out  as  an 
addition  to  tha  initial  rain  data-reduction  task  because  the  drop-size 
distributions  were  thought  to  be  the  key  to  the  problem  of  the  large 
observed  variations  in  backscatter  for  a constant  rain  rate.  The 
Illinois  Water  Survey  Group,  which  were  originally  contracted  to  analyze 
the  drop-size  tapes,  were  unable  to  do  sc  due  to  the  generally  poor 
quality  of  the  tapes.  Because  of  the  importance  of  the  drop-size  dis- 
tributions to  the  backscatter  question,  EES  proposed  to  attempt  the 
analysis  of  the  tapes  as  an  additional  task  to  the  initial  contract. 

The  rain  drop-size  tapes  were  thus  obtained  from  the  Illinois  Water 
Survey  and  processed  to  recover  as  such  of  the  drop-size  data  as  possible. 

The  drop-size  spectrometers  utilize  a piezoelectric  crystal  to 
convert  the  impact  of  raindrops  on  the  spectrometer  heads  to  voltage 
pulses  which  are  then  record-d  on  a cassette  recorder.  The  spectrometers 
are  calibrated  by  dropping  plastic  beads  of  various  diameters  onto  the 
sensor  heads  from  a tower  and  .ecording  the  outputs.  By  this  calibration 
technique,  the  Illinois  tfa'vr  Survey  was  able  to  determine  that  the 
relationship  between  drop  diameter  ».nd  the  peak  of  the  corresponding 
voltage  pulse  was  approximately  a .linear  function  for  these  spectrometers. 

In  order  to  account  for  changes  in  recorder  gain  and  offset  over  a period 
of  time,  at  the  onset  of  each  rain  storm,  a calibration  signal  was  auto- 
matically recorded  on  the  tape  consisting  of  a 100-Hz  square  wave  with  an 
amplitude  equal  to  the  height  of  a voltage  pulse  corresponding  to  the  impact 
of  a 6-mc  rain  drop.  In  addition,  time  code  uas  recorded  on  a second 
record  channel,  but  unfortunately  this  feature  did  not  function  most  of 
the  time. 
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Figure  11  gives  a block  diagram  of  the  equipment  set-up  used  to 
analyze  the  drop-size  data.  The  output  of  the  data  recorder  was  passed 
through  a variable  gain  and  offset  amplifier  to  a peak  detector.  The 
peak  detector  was  required  because  the  output  of  the  piezoelectric  sensor 
exhibited  an  under-damped  sin(x)/x  type  of  response  following  the  impact 
of  a raindrop,  whereas  only  the  peak  of  the  signal  is  of  interest.  The 
peak  detector  was  triggered  by  a Tektronix  545A  oscilloscope  which 
served  as  a variable  threshold  trigger  generator.  When  a rain  drop  im- 
pacted, the  main  sweep  of  the  oscilloscope  was  triggered,  which  in  turn 
triggered  a pulse  generator  that  enabled  the  peak  detector  for  250  psec 
following  the  start  of  the  Tektronix  scope  main-sweep.  The  Fabritek 
computer  was  used  in  the  pulse-height  analyzer  mode  (see  discussion  of 
the  radar  backscatter  amplitude  data  reduction)  to  generate  a density 
function  consisting  of  the  number  of  occurences  versus  drop-size  (to  the 
nearest  0.1  nsa)  during  a given  time  interval  (usually  30  seconds). 

This  distribution  was  then  utilized  to  calculate  both  the  rain  rate 
and  a normalized  distribution  function  consisting  of  the  number  of  drops, 

of  a given  diameter,  that  were  contained  in  a cylinder  above  the  sensor 

3 

head  of  the  drop-size  spectrometer  with  a volume  of  1 ta  . The  rain 
rate  can  be  expressed  as  follows: 

R(mn/lir)  * 3600  H/T, 

where  H is  the  height  of  accumulated  water  over  the  sensor  head  in  mo 
and  T is  the  time  in  seconds. 

H is  equivalent  l > the  v ilume  of  th<"  water  divided  by  the  sr  sor 
2 

head  area  (3848  na  ).  The  volume  is  given  by: 


n 


i * 1 


where:  is  the  number  of  drops  of  a given  diameter  (+  0.05  mm) 

which  strike  the  sensor  head,  and 

D is  the  maximum  drop  diameter, 
n r 

Combining  the  above  factors  yields: 
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The  volunetric  distribution  (nuaber  of  drops  in  a cylindrical  coition 
above  the  sensor  head  with  voluae  of  one  cubic  Deter)  is  related  to  the 
measured  distribution  over  some  tine  Interval  by: 

Vih 

N (vol.)  * ■ T-  — 

where:  is  the  nuaber  of  drops  striking  the  sensor  in  tine  T with 

diameter 

h is  the  height  of  a cylindrical  colunn  with  the  sane  diameter  as 

3 

the  sensor  head  and  a voluae  of  1 & (260  meters  for  these  sensors) , 

and  v is  the  terminal  velocity  of  a drop  diameter  D^. 

The  terminal  velocities  used  to  calculate  the  volunetric  distributions 
were  taken  fron  Gunn  and  Kinzer  [6]. 

A Focal  program  was  written  for  the  PDP-8/F  conputer  to  perform  the 
above  calculations,  and  to  allow  graphing  of  the  rain  rates  and  voluae trie 
drop-size  distributions  on  an  X— Y plotter.  The  drop— size  distributions  were 
determined  -s  closely  as  possible  for  the  sane  time  intervals  as  were 
analyzed  radar  backscatter  data  in  order  to  determine  any  relationships. 

A nuaber  of  problens  were  encountered  in  the  data  reduction.  The 
tape  recorders  used  originally  to  record  the  data  were  not  of  adequate 
quality;  due  to  poor  frequency  response,  the  calibrations  cn  the  ta >es 
were  often  so  distorted  as  to  be  unusable.  Instead,  the  distributions 
were  calibrated  by  comparing  the  calculated  rain  rates  from  the  drop-size 
distributions  tc  the  tipping  bucket  records.  Although  this  was  a tedious 
process,  it  appeared  to  work  quite  well,  as  the  match-up  of  rain  rates  be- 
tween the  spectrometers  and  the  tipping  buckets  in  general  was  quite  good. 
The  lack  of  tine  code  on  aost  of  the  tapes  further  complicated  the  process, 
but  once  again  the  tipping  bucket  records  were  used  to  determine  times  for 
the  drop-size  spectrometer  tapes.  In  addition,  high  noise  levels  on  the 
tapes  prevented  the  detection  of  some  of  the  small  drops. 
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However,  a number  of  distributions  were  obtained  which  correlated  well 
in  rates  with  the  tipping  buckets  and  appeared  reasonable  so  as  to  allow 
sene  general  determinations  of  the  relationship  between  drop-size  and 
radar  backscatter  characreristics  to  be  sade.  These  findings  will  be 
discussed  in  a later  section. 

3.  Sources  of  Error  in  the  Data 

The  analyzed  backscatter  data  show  large  variability  in  both  amplitude 
and  spectral  parameters  when  plotted  as  a function  of  rainfall  rate. 

While  this  probably  indicates  that  rainfall  race  is  not  an  adequate  and 
sufficient  parameter  for  backscatter  characterization,  there  were  some 
problems  in  the  experiment  and  in  the  quality  of  the  data  on  the  tapes 
which  could  have  affected  the  accuracy  of  '.he  analyzed  data. 

There  were  two  main  problems  with  the  way  in  which  the  experiments 
were  structured.  (1)  The  rain  instrumentation  was  placed  near  the  radar 
targets,  while  data  on  the  backscatter,  of  necessity,  had  to  be  sampled 
at  radar  cells  which  did  not  include  the  radar  targets.  Since  the  rains 
seen  during  the  field  tests  were  often  nonuni  forts,  the  rain  rates  and 
drop-sizes  as  measured  by  the  rain  instruments  may  not  in  general  be  the 
same  as  those  in  the  radar  cells  being  sampled.  (2)  The  signal  back- 
scattered  from  rain  is  attenuated  when  propagating  through  a rain-filled 
volume.  Because  of  the  way  the  experiment  was  instrumented,  there  was  no 
direct  way  to  measure  this  attenuation  for  the  data  recorded  on  tape, 
so  that  an  attenuation  function  had  to  be  assumed  from  theoretical  cal- 
culations. This  fact,  coupled  with  the  problem  that  many  of  the  measurements 
were  made  at  ranges  which  were  in  the  near  field  of  the  70-GHz  and  95-CHz 
antennas,  could  have  contributed  further  to  errors  in  the  backscatter. 

Several  problems  in  the  recorded  data  also  could  have  caused  errors 
in  the  analyzed  data.  Several  of  the  tapes  had  no  calibrations  recorded 
on  them,  so  that  calibrations  had  to  be  used  from  tapes  recorded  on 
different  days.  Obviously,  any  changes  in  radar  parameters  such  as  trans- 
mitted power  or  reciever  drift  could  result  in  calibration  errors. 
Furthermore,  the  35-GKz  radar  had  an  extremely  nonlinear  transfer-function, 
in  terms  of  d3,  which  could  have  increased  the  errors  in  the  data  analysis. 
Xorcally  the  RMS  error  in  processing  the  recorded  data  is  about  + 0.5  dB, 
but  could  have  been  higher  in  this  case.  Also  quite  large  60-Hz  and  400-Hz 
components  were  present  in  much  of  the  data,  requiring  the  use  of  notch 
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filters.  This  was  particularly  a problea  in  attempts  to  jnaJ  the 
spectral  shape  of  the  data,  since  the  notch  filters  distorted  the 
spectrua.  As  explained  previously,  this  problen  was  partially  solved 
by  the  use  of  an  algorithn  in  the  PDP-8/F  cocputer  which  oultiplied  the 
filtered  spectrua  by  the  filter  inverse  functioa  after  deletion  of  the 
60-Hz  or  4 00- Hz  spikes. 

It  can  therefore  be  seen  that  the  spread  in  the  data  can  be  explained 
to  soce  extent  by  the  problems  encountered  in  the  data  analysis.  However, 
the  data,  in  spite  of  the  possible  errors,  does  sees  to  be  consistent  enough 
to  allow  conclusions  as  to  the  aechanisns  involved  in  the  scattering,  and 
to  perait  determination  of  worst-case  conditions  for  backs cat ter  amplitude 
and  spectral  width. 

B.  Suacary  of  Results 

1.  Interpretation  of  Data, 

Periodic  waves,  such  as  radar  signals,  incident  upon  a nateriai  body 
can  be  described  in  terns  of  absorbed  energy  and  scattered  energy.  The 
scattered  wave  is  obtained  by  solving  Maxwell's  equations  at  the  boundary 
of  the  body,  where  the  geoaetrical  description  cf  the  body  thus  dictates 
the  coordinate  systen  and  the  natheoatical  functions  required  for  a 
solution.  For  the  case  of  cylindrical  and  spherical  objects,  the  Bessel 
functions  and  Legendre  functions  are  required  to  describe  the  scattered 
waves.  The  variable  paraaeter  in  these  functions  is  the  normalized 
particle  aiaaeter  in  units  of  wavelength. 

When  the  argunent  of  the  Bessel  function  is  snail,  i.e. , if  the 
radius  of  curvature  is  auch  smaller  than  the  wavelength  of  the  incident 
signal,  then  only  the  first  ter®  in  the  series  representation  cf  the 
Bessel  function  is  usee.  The  resulting  solutions  of  Maxwell's  equations 
take  the  sane  fora  as  used  by  Rayleigh  in  describing  scattering  of  light 
fros  snail  particles.  An  effective  cross-section  Is  defined  in  tens  of 
the  ratio  of  reflected  power  to  incident  power  and  the  equation  takes 
the  fora 
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radar  cross-section  and  attenuation  to  increase  as  rain  rate  increases. 

The  exact  values  vary  from  one  investigator  to  another,  depending  on  the 
models  employed. 

Rain  rate  and  drop-size  distributions  have  a general  relationship  if 
measured  over  long  time  intervals  and  averaged  over  many  minutes.  However, 
it  has  been  well  documented  that  a wide  variation  in  the  total  number  of 
drops  and  in  drop-size  distribution  will  occur  over  intervals  of  a few 
seconds  to  a minute  at  a given  rain  rate  [8]. 

2.  Attenuation 

Attenuation  due  to  uncondensed  gases,  principally  water  vapor  and 
oxygen,  ware  not  considered  important  during  these  tests  since  the  range 
was  so  short.  At  70  GHz,  the  highest  attenuation  for  any  of  the  frequencies 
used,  the  attenuation  due  to  uncondensed  gases  is  only  a few  tenths  of  a 
dE  at  worst. 

However,  attenuation  caused  by  the  water  droplets  themselves  during 
rain  storms  is  an  important  parameter.  The  data  analyzed  under  this  pro- 
ject were  obtained  from  range-sampled  signals  and  thus  the  signal  strength 
is  affected  by  the  amount  of  path  attenuation.  During  analysis,  it  was 
not  possible  to  obtain  values  of  attenuation  since  range  samples  were  not 
recorded  under  constantly  definable  conditions.  A search  was  made  of  the 
literature  to  obtain  the  most  reasonable  average  value  of  attenuation 
caused  by  rain  at  the  radar  frequencies  of  interest. 

Calculations  made  by  J.  de  Bettencourt  at  Raytheon,  [9]  by  Mueller 
and  Sims  at  Illinois  State  Water  Survey  [10],  and  by  Lin  and  Ishimaru  at 
the  University  of  Washington  [11]  appear  to  be  the  most  reasonable  as 
applied  to  this  project.  Table  2 contains  the  attenuation  due  to  rain, 
as  a function  of  rain  rate  and  radar  frequency,  which  have  been  used 
during  the  analysis  efforts  on  this  project. 

3.  Average  Return  Amplitude 

The  amount  of  backscatter  per  unit  volume  of  radar  resolution  has 
been  calculated  and  is  presented  in  Figures  12  and  13  for  two  polariza- 
tions. These  data  could  be  compared  with  the  many  values  presented  in  the 
literature  and  suggested  by  Richard  and  Kammerer  [3],  and  L.  D.  Strom  [12], 
Basic  experimental  data  used  for  the  results  presented  in  this  report  were 
collected  by  the  same  operators  as  those  collecting  dat.  for  Richard  and 
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TABLE  2 


Attenuation  Due  to  Rain 


f = 
A = 

a = 

10  GHz 
3.2  cm 

0.00919R1*16 

35  GHz 
0.86  cm 

0.273R0*985 

70  GHz 
0.43  cm 

0.63<R°-868 

95  GHz 
0.32  cm 

1.6R0*64 

R 

(mm/hr) 

a (dB/Km) 

0 

0 

0 

0 

0 

5 

0.509 

1.33 

2.56 

4.48 

10 

0.133 

2.64 

4.68 

6.984 

15 

0.213 

3.93 

6.65 

9.05 

20 

0.297 

5.22 

8.54 

10.88 

25 

0.385 

6.50 

10.36 

12.55 

30 

0.475 

7.78 

12.14 

14.11 

35 

0.568 

9.06 

13.88 

15.57 

40 

0.663 

10.33 

15..8 

16.96 

45 

0.760 

11.60 

17.26 

18.29 

50 

0.859 

12.87 

18.91 

19.56 

55 

0.960 

14.14 

20.55 

20.79 

60 

1.062 

15.40 

22.16 

21.99 

65 

1.165 

16.67 

23.75 

23.14 

70 

1.270 

17.93 

25.32 

24.27 

75 

1.375 

19.19 

26.89 

25.36 

80 

1.482 

20.45 

28.44 

26.43 

85 

1.590 

21.71 

29.98 

27.48 

90 

1.700 

22.97 

31.50 

28.50 

95 

1.809 

24,22 

33.02 

29.50 

100 

1.920 

25.48 

34.52 

30.49 

Mueller-Sims  1969  * 
Florida  Rain 

J,  de  Bettencourt 
1973  Raytheon 
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Rain  Rate  (ran/hr) 


Figure  12.  Least  mean  square  fit  to  radar  backscatcer  data 
W polarization. 


Kararaerer,  as  was  discussed  in  a previous  section. 

Excellent  agreement  between  data  presented  here  and  the  previously 
reported  data,  both  experimental  and  theoretical,  exists  at  X-band.  This 
is  gratifying,  since  one  of  the  primary  reasons  for  having  an  X-band 
measurement  system  was  to  establish  the  data-base  reference  and  validate 
the  entire  experiment.  Good  agreement  between  theory  and  experiment  was 
expected  since  drop-size  distributions  show  the  maximum  drop  diameter  to 
be  less  than  8 me;  thus,  at  X-band  (A  = 22  era)  the  scattering  mechanism 
is  Rayleigh. 

On  the  average,  the  backscatter  obtained  at  radar  frequencies  of 
9.374  and  35  GHz  are  parallel  to  each  other,  whic..  indicates  that  Rayleigh 
Scattering  is  the  primary  mechanism  at  both  those  frequencies  However, 
it  is  noted  that  the  backscatter,  as  a function  of  rain  rate,  o.  :omes 
flatter  at  70  and  95  GHz  when  compared  to  the  lower  frequencies.  This 
might  be  expected  since  the  radar  wavelength  is  comparable  to  rain  drop 
size  and  the  mechanism  is  Mie  Scattering.  This  same  trend  occurred  for 
both  vertical  polarization  and  right  circular  polarization. 

In  general,  circular  polarization  produced  a lower  backscatter  than 
did  linear  polarization.  Differences  obtained  during  this  investigation 
are  smaller  than  reported  previously.  A fairly  consistent  10  to  15  dB 
difference  is  reported  here,  whereas  some  investigators  have  reported 
circular  returns  20  dB  lower  than  linearly  polarized  returns. 

An  Interesting  result  appeared  in  the  averaged  calculated  backscatter 
produced  on  this  phase  of  the  program  but  not  observed  by  other  inves- 
tigators. The  backscatter  at  95  GHz  was  lower  than  at  70  GHz  for  both 
linear  and  circular  polarization.  It  is  postualted  that  some  forward 
scattering  occurred  at  95  GHz.  This  concept  could  be  validated  by  measure- 
ments at  140  GHz,  for  which  the  transmitted  signal  would  have  a wave- 
length comparable  to  the  rain  drop  size.  If  a resonant  phenomenon  is 
creating  forward-scattering,  the  backscatter  at  140  GHz  should  be  even 
lower  than  at  95  GHz. 

4.  Amplitude  Fluctuations 

As  discussed  in  a previous  section,  the  radar  signal  was  selected  by 
a range  gate  sampler  at  selected  ranges.  The  signal  recorded  on  magnetic 
tape  thus  represents  the  return  from  a volume  defined  by  the  antenna  beam 
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width  and  the  transmitted  pulse  length.  Considering  all  the  various  para- 
meters of  the  four  radar  systems  and  the  various  ranges,  the  volumes  spanned 
from  5 to  10,000  cubic  meters. 

From  previously  reported  date  [8]  [10],  the  number  of  drops  in  a cubic 
meter  varies  as  the  rain  rate  increases*  In  general,  for  diameters  around 
0.3  mm,  the  concentration  is  100  to  200  drops  per  cubic  meter,  almost 
independently  of  rain  rate.  However,  as  the  rain  rate  increases,  larger 
sized  drops  occur  with  maximum  diameters  up  to  8 ran.  As  the  rain  rate 
increases  from  4 ran/hr  to  16  ran/hr,  on  the  average  the  number  of  drops, 
greater  than  0.5  ran  in  diameter  will  increase  from  a few  drops  to  several 
hundred,  and  the  mean  diameter  will  increase  from  around  1.4  to  2.4  tan. 

From  the  system  parameters,  the  test  area  geometry,  and  previous 
rain  data,  it  was  anticipated  that  during  light  to  moderate  rain  rates 
there  should  be  a large  number  of  scatterers  and  drops  of  sufficient 
size  to  produce  measurable  backscatter  radar  signals.  This  was  indeed 
observed  during  the  field  data  collection  period.  Analysis  of  the 
recorded  data  shows  that  the  amplitude  of  the  return  signal  fluctuates 
in  a log-noraal  manner.  In  addition,  analysis  reveal  that  this  log- 
normal characteristic  exists  at  all  rain  rates  from  1 ro/hr  to  100  cxi/hr 
and  at  all  the  operating  frequencies  from  9.375  Hz  to  95  GHz.  Examples 
of  probability  densities  and  cumulative  distributions  for  the  backscatter 
from  rain  are  shown  in  Figures  14,  15,  16,  and  1.'. 

Calculation  of  the  standard  deviation  and  variance,  from  the  log- 
normal distribution  shows  that  for  any  of  the  radar  systems,  the  standard 
deviation  is  independent  of  rain  rate.  This  implies  that  frcn  lipht 
rains  to  extremely  heavy  rains  that  there  were  a sufficient  number  of 
rain  drops  having  a random  location  and  distribution  and  without  phase 
coherence  within  the  rain  cell  resolution,  sj  that  the  total  s-gnal  re- 
flected back  to  the  radar  was  normally  distributed.  The  probability  plots 
obtained  from  right-circular  and  vertical  polarization  are  both  log- 
normal. The  mean  values  for  circular  polarizations  were,  of  course,  con- 
sistently lower  than  mean  values  for  linear  polarizations. 

Standard  deviations  obtained  from  vertical  polarizations  are  shewn 
in  Figures  18  through  21  and  for  circular  polarization  in  Figures  21 
through  25. 
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Figure  16.  Probability  density  end  cumulative  distribution*  for  backscatter  from  rain; 
70  GHz,  5 mm/hr  rain  rate,  W polarization. 
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i^uro  20.  Rain  backscatter  amplitude  standard  deviation  versus  rain 
rate  for  70  Qlz  frequency,  W polarization. 
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i'j ;;i«re  21.  Rain  backscatter  amplitude  standard  deviation  versus  rain 
rate  lor  95  GHz  frequency,  W polarization. 


5.  Power  Density  Spectrum 

All  the  radars  employed  in  these  tests  were  incoherent  systems.  Thus, 
it  was  not  possible  to  obtain  the  coherent  Doppler  frequency  spectrum  of 
the  return  signals  from  the  rain  clutter.  However,  it  is  possible  to 
determine  the  frequency  components  without  coherent  phase  relationship  in 
the  power  spectrum.  The  equipment  and  data-processing  techniques  for 
obtaining  the  spectra  are  presented  in  a previous  section  of  this  report. 

For  the  purpose  of  comparing  data  from  this  project  to  that  obtained 
by  previous  investigators,  it  is  important  to  note  differences  in  the 
measurement  radars  employed.  Three  of  the  measurement  radars  used  in  this 
investigation  employed  log-if  amplifiers;  the  35-GHz  (APQrl37)  system  uses 
a lin-1cg  if  amplifier.  Previous  investigators  have  used  linear  receiver 
systems.  The  analysis  technique  of  using  fast  Fourier  transforms  to 
convert  time-domain  signals  to  frequency-domain  signals  requires  that 
the  signal  be  properly  formated  prior  to  computer  processing. 

Since  many  current  radar  systems  are  designed  with  logarithmic 
receivers,  it  was  decided  to  process  the  data  in  the  same  form  as  they  were 
recorded.  This  will  allow  system  designers  to  use  the  final  results 
directly.  Thus,  if  a Doppler  processor  is  used  with  log  receivers,  then 
the  spectral  width  expected  from  rain  backscStter  will  be  as  shown  in 
Figures  26  through  25  for  vertical  polarization,  and  in  Figures  30  through 
31  for  circular  polarization. 

The  frequency  bandwidth  increases  as  rain  rate  increases,  which  is 

reasonable  since  having  more  particles  allows  more  reflections  and 

thus  spreads  out  the  spectrum  in  the  frequency  domain.  It  is  also  noted 

that  as  the  radar  frequency  was  increased,  the  spectrum  width  increased, 

but  the  amplitude  of  the  spectrum  decreased.  This  is  consistant  with  the 

decrease  in  the  standard  deviation  of  the  return  when  comparing  X-band 

to  95  GHz  signals.  When  observing  the  spectrum  at  a high  rain  rate 

(greater  than  35  mm/hr)  for  the  70  and  95  GHz  radars,  the  return  appears 

almost  as  colored  noise.  It  should  be  noted  that  the  roll-off  of  the 

spectrum  with  frequency  is  not  Gaussian,  but  rather  is  better  fitted  by 

a 1/(1  + ( f / f ) response  curve  at  10  GHz  (which  is  similar  to  the 

spectrum  shape  of  the  return  from  land  clutter),  and  is  more  closely 

2 

represented  by  a 1/(1  + f/f()  ) response  at  35  GHz  and  higher.  Thus, 
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VV  polarization 


formalized  frequency  spectrum  of  rain  return  at  95  GHz 


the  assumption  of  a Gaussian  shaped  bandpass  for  rain  clutter  spectra 
for  a radar  system  design  could  lead  to  serious  errors  in  the  calculated 
spectral  bandwidth. 

On  a number  of  spectra,  a second  distribution  appears  removed  from 
the  low-frequency  continous  distribution.  Application  of  the  standard 
Doppler  equation, 

f - 102  x velocity 
wavelength 

reveals  that  the  center  frequency  was  produced  hy  the  movement  of  the  total 
rain  cell  in  a radial  direction  from  the  radars.  For  every  case  tested, 
the  field  logs  indicate  a wind  speed  and  direction  which  would  produce  a 
Doppler  frequency  within  10a  of  that  calculated.  This  accuracy  was  much 
higher  than  expected;  however,  it  has  been  verified  on  enough  occasions 
that  one  is  led  to- believe  that  the  data  reduction  process  is  accurate. 

(For  an  example  of  this  type  of  spectrum,  see  Appendix  b. ) 

6.  Correlation  Functions 

Of  concern  to  the  designer  of  radars  used  in  a space  scanning  node, 
is  the  time  required  to  obtain  independent  samples  of  radar  returns. 

Tiie  auto-correlation  function  provides  this  basic  information  so  that 
the  radar  backscatter  data  were  processed  to  obtain  the  linear  correlation 
functions.  The  results  of  this  investigation  are  shown  in  Figure  32. 

Each  point  on  this  figure  is  the  result  of  8 data  runs  each  lasting 
approximately  0.7  seconds.  It  is  noticed  that  in  general  as  the  rain 
rate  increases  and/or  the  radar  frequency  is  increased,  the  decorrelation 
time  decreases.  This  is  consistant  with  the  observation  that  the  power 
spectrum  distribution  becomes  more  noise  like  as  the  rain  rate  and/or  the 
radar  frequency  is  increased. 

Values  obtained  from  this  investigation  are  consistant  with  those 
reported  by  Goldhirsh  and  Katz  (13J.  They  reported  that  on  the  gross 
average  one  can  expect  decorrelation  times  as  follows:  S-band,  30  milliseconds 
X-band  10  milliseconds;  and  K-band,  3 milliseconds.  Examples  of  specific 
auto-correlation  functions  are  shown  in  Figures  33  and  34. 

Cross  correlations  functions  were  obtained  between  the  signals  from 
the  four  radar  systems  when  operating  at  exactly  the  same  time.  Regardless 
of  rain  rate  or  transmitted  polarization,  the  signals  between  the  radars 
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Decorrelation  time  vereus  tain  rate  and  frequency 
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Figure  33.  Normalized  autocorrelation  function  for  rain  backscatter 
9.375  GHz  frequency,  9 k/ hr  rain  rate,  W polarization. 
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Figure  34.  Normalized  autocorrelation  function  for  rain  backscatter 
95  GHz  frequency,  9 na/hr  rain  rate,  W polarization. 
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were  decorrelated  z'c  all  elves.  Since  the  cross-correlation  coefficients 
were  always  zero , no  graphs  are  presented. 


The  data  froe  rain  backseat ter  as  a function  of  rain  rate,  as  shown 
in  Figures  12  and  13,  were  derived  from  rain  tipping  buckets.  These 
devices  are  aeasuresent  units  which  average  rain  rate  over  a relatively 
long  interval  of  time  as  coopared  to  the  observation  period  of  a radar 
system.  The  buckets  used  during  this  experiment  have  a switch  closure 
every  .01  inch  (0.254  eh)  of  water  entering  the  bucket.  Thus  if  the  rain 
rate  is  2.54  ma/hr,  there  will  be  10  closures  every  hour,  or  one  closure 
every  6 minutes  while  the  radar  data  were  evaluated  for  30-second  intervals. 
If  good  rain  rate  sampling  were  performed,  then  there  should  be  two 
samples  of  the  rain  rate  during  the  30  second  radar  sample  period.  With 
the  tipping  buckets  employed,  this  would  require  a minimus  rain  rate  of 
61  om/hr.  Obviously  there  is  a conflict  and  the  result  is  that  a "gross- 
averaging"  effect  was  involved  whenever  rain  rate  was  related  to  any 
particular  radar  observation. 

An  example  can  be  made  by  looking  at  Figure  35.  There  was  a total 
lapsed  time  of  2 minutes,  in  which  the  radar  returns  were  evaluated  4 times 
at  30-second  intervals.  The  rain  drop-size  spectrometer  data  were  also 
evaluated  for  four  20-second  intervals,  and  the  rain  rates  were  calcu- 
lated from  the  drop-size  distributions.  The  tipping  buckets  indicated 
an  average  rain  rate  of  2.7  m/hr,  whereas  the  spectrometers  show  a dynamic 
change  in  rate  varying  from  4.5  m/hr  to  9.3  m/hr,  with  an  average  of 
7.25  m/hr  over  this  short  2-ainute  interval.  The  difference  in  average 
rain  rate  as  determined  by  the  tipping  buckets  compared  to  the  spectr- 
meters  can  be  accounted  for  by  the  fact  that  the  tipping  bucket,  at  this 
rain  rate,  was  averaging  over  5 to  6 minute  intervals.  (This  was  the  actual 
time  between  tips.) 

In  addition  to  rapidly  changing  rain  rates,  dynamic  changes  in 
drop-size  distribution  were  also  common.  Since  drop-size  plays  a major 
role  in  the  magnitude  of  radar  backscatter,  this  Is  a source  of  back- 
scatter  variance,  as  discussed  earlier. 

Due  to  a variety  of  equipment  problems,  not  all  radar  data  were 
recorded  on  magnetic  tapes,  and  all  that  were  recorded  could  not  be 
analyzed.  In  addition,  not  all  the  rain  spectrometers  were  functional 
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It  has  been  observed  that  at  low  rain  rates,  the  tipping  buckets 
employed  did  not  respond  fast  enough  to  present  a real  time  measure  of 
the  rain  rate.  This  affects  not  only  the  value  of  backscatter  but  also 
the  path  attenuation.  Thus  it  is  not  unexpected  that  a wide  spread  in 
data  points  resulted  at  any  given  rain  rate. 


III.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  summarizes  results  of  a thorough  analysis  of  data  obtained 
during  a series  of  measurements  of  radar  backscatter  from  rain  at  fre- 
quencies of  9.375,  35,  70,  and  95  GHz.  The  geometry  of  the  experiment  and* 
the  equipment  used  were  chosen  with  the  intent  to  provide  data  that  would 
be  useful  to  an  equipment  designer  in  specifying  an  optimum  operating 
frequency  for  a given  mission.  It  is  believed  that  this  goal  has  been 
achieved;  however,  a number  of  experimental  difficulties  were  encountered 
which  complicated  the  analysis,  and  which  suggest  that  for  future  experiments 
in  this  area  special  care  should  be  given  to  the  choice  and  use  of  support- 
ing instrumentation.  A particularly  important  requirement  exists  for  more 
extensive  use  of  meterological  instrumentation  which  is  reliable  and  whi'ih 
can  be  readily  calibrated.  Not  only  is  extensive  rain  rate  and  drop-size 
distribution  measuring  equipment  required,  but  also  other  classes  of 
instrumentation  for  measuring  such  parameters  as  wind  speed,  wind  direction, 
and  air  temperature. 

The  emphasis  of  the  measurements  and  subsequent  analysis  was  directed 
toward  investigating  the  amplitude  and  noncoherent  frequency  characteristics 
of  the  radar  return.  Limited  checks  of  the  attenuation  characteristics 
of  the  path  were  attempted;  however,  since  the  measurements  were  not 
specifically  designed  to  allow  accurate  determination  of  loss,  only  general 
agreement  with  published  data  was  obtained.  Because  of  limitations  of  the 
loss  calculations  based  on  the  data  reported  here,  detailed  data  from 
the  literature  have  been  used  to  determine  the  approprate  attenuation 
values  as  described  below.  It  is  recoanended  that  future  experiments  of 
this  type  be  designed  to  address  the  question  of  path  loss  directly  as 
one  of  the  main  objectives. 

Rain  rates  qccuring  during  the  test  program  varied  frc  less  than  1 
mm/hr  to  more  than  90  ra/hr.  The  average  backscatter  for  linear  polari- 
zation was  observed  to  increase  linearly  (on  a log-log  plot) , with  in- 
creasing rain  rate  at  all  four  frequencies.  When  circular  polarization 
was  used,  the  returns  were  consistently  10  to  15  dB  less  than  when  using 
linear  polarization.  The  values  obtained  at  X-band  and  Ka-band  are 
accurately  represented  by  theoretical  calculations  using  Rayleigh 
scattering  and  drop-size  distributions  according  to  Laws  and  Parsons  [14]. 
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The  backscatter  values  obtained  at  70  and  95  GHz  have  a smaller  slope  than 
at  S.375  and  35  GHz.  This  is  generally  typical  of  the  change  expected  when 
scattering  is  predominately  in  the  Mie  region.  Since  Mie  scattering  is 
extremely  sensitive  to  the  exact  nature  of  the  drop-size  distribution , a 
number  of  actual  distributions  were  checked  and  found  to  confirm  both 
the  general  change  in  slope  observed  and  the  significant  scatter  of  the 
average  cross-section  observed  between  runs.  One  anomaly  noted  here  was 
that  the  backscatter  at  95  GHz  on  the  average  was  less  than  that  at  70 
GHz.  It  is  suggested  that  this  may  result  from  a reasonance  between 
drop-size  and  radar  wavelength,  thereby  supporting  the  concept  of  en- 
hanced forward  scattering  at  these  wavelengths.  This  hypothesis  could 
not  be  checked  here  as  the  experiment  did  not  provide  the  appropriate 
comparison  data. 

The  rain  storms  observed  during  these  measurements  varied  strongly 
in  spatial  and  temporal  nature,  but  were  generally  of  relatively  short 
duration.  These  characteristics  were  major  contributors  to  the  extreme 
variability  of  the  drop-size  distributions  observed  and  thus,  consequently, 
the  variability  of  the  backscatter.  It  is  believed  that  the  extreme 
variability  of  the  storms  observed  made  them  represent  essentially  a 
"worse  case"  condition  both  from  an  experimental  point  of  view  and  a 
radar  systems  design  viewpoint. 

Amplitude  distributions  of  the  backscatter  were  observed  to  be 
approximately  log-normal  at  all  frequencies  and  for  rain  rates  between 
2 and  60  mrn/hr.  For  rain  rates  less  than  1 an /hr  and  greater  than 
60  nm/hr,  the  amplitude  distributions  were  not  consistently  log-normal 
but  were  generally  skewed.  This  may  have  been  caused  by  having  a low 
signal-to-noise  ration,  or,  at  the  higher  return  levels,  by  system 
saturation;  however,  it  is  also  possible  that  the  skewness  is  a real 
phenomenon  due  to  the  observed  behavior  of  the  drop-size  distributions 
at  the  extremes.  In  addition,  the  amplitude  variance  of  the  returns 
appeared  to  be  independent,  on  the  average,  of  rain  rate,  and  to  generally 
decrease  with  increasing  frequency.  The  variance  was  generally  less  for 
circular  polarization  than  for  the  linear  polarization. 

The  noncoherent  frequency  spectral  width  was  observed  to  increase 
with  rain  rate  at  each  test  frequency,  and  also  to  increase  with  frequency 
at  a given  rain  rate.  The  spectrum  amplitude  was  also  observed  to  decrease 
with  increasing  frequency,  approaching  noiselike  characteristics  at  95  GHz. 
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The  roll-off  rate  at  the  high-frequency  end  of  the  spectrum  was  observed 
to  be  much  slower  than  Gaussian,  being  best  described  as  quadratic  in 
nature.  The  spectra  at  X-band  were  generally  well  described  by  a cubic 
function.  No  significant  correlation  (as  obtained  from  the  cross- 
correlation function)  was  discernible  between  the  returns  from  the  four 
radars  when  observed  at  the  same  time.  The  decorrelation  time  (as 
obtained  from  the  auto-correlation  function)  became  less  with  increasing 
rain  rate  and  with  increasing  radar  frequency.  Values  of  the  de- 
correlation time  ranged  from  14  milliseconds  at  5 m/hr  r3in  at  X-band 
to  1.4  milliseconds  at  100  mm/hr  rain  at  95  GHz. 

It  is  concluded  that,  while  much  new  data  are  available  to  the 
designer  as  a result  of  this  work,  a number  of  important  questions  remain 
unanswered  as  to  the  impact  of  rain  on  system  performance.  Perhaps  the 
hardest  question  to  settle  satisfactorily  is  still  that  of  the  compara- 
tive merits  of  a millimeter  radar  with  other  approaches  (i.e.  microwave 
or  laser  radar)  in  moderate  rain  and/or  fog,  for  precision,  short-range 
applications  in  situations  where  space  or  geometry  limit  the  size  of 
antenna.  Examples  of  such  situations  include  helicopter  systems  and 
most  point-defense  applications.  Data  reported  here  suggest  that  95 
GH2  is  very  competitive  in  such  applications.  Experiments  should  be 
conducted  which  make  the  necessary  comparisons  between  the  competing 
approaches  in  rain,  while  giving  adequate  attention  to  the  required 
meteorological  conditions  to  resolve  this  question. 
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APPENDIX  A 


Plotted  Average  Backscatter  Data 

This  appendix  contains  the  plotted  values  of  average  radar  back- 
scatter data  versus  rain  rate  which  were  used  to  calculate  the  least- 
squares  functions  presented  in  Figures  12  and  13.  Each  point  represents 
* approximately  30  seconds  of  data  which  were  measured  at  one  of  four 

different  ranges.  The  values  have  been  normal iced  to  remove  dependence 

on  range  and  systen  parameters , and  thus  have  units  of  radar  cross- 

2 3 

section  per  unit  volume  (a  /q  ).  As  has  been  previously  noted,  the 
data  indicate  that  variations  of  greater  than  10  dB  in  the  magnitude 
of  the  backscatter  can  occur  for  a constant  rain  rate.  As  discussed 
in  Section  III-B-7,  changes  in  the  drop-size  distribution  constitute 
one  important  cause  for  such  large  variations  in  backscatter. 
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APPENDIX  B 

Spectral  Distributions  for  Selected  Data  Runs 

This  appendix  contains  selected  examples  of  spectral  plots  of  the 
fluctuations  of  received  backscatter  from  rain  for  frequencies  of 
9.375,  35,  70,  and  95  GHz,  rain  rates  of  5 to  100  oo/hr,  and  vertical 
and  circular  polarizations.  Since  all  the  radars  utilized  logarithmic 
receivers,  the  frequency  spectra  of  their  output  voltages  represent 
fluctuations  about  the  mean  values  and  thus  have  amplitudes  calibrated 
in  dB. 

These  examples  illustrate  the  general  trends  observed  for  all  the 
data,  that  is  for  increasing  frequency  and  rain-rate,  the  peak  spectral 
amplitude  of  the  fluctuations  decreases  and  the  spectral  width  increases. 

In  the  limit,  for  95  GHz  frequency  and  100  ao/hr  rain-rate,  the  spectrum 
appears  very  noise-like  with  less  than  1 dB  of  fluctuation  and  a very 
wide  spectral  width.  The  spectra  also  appear  polarization-dependent 
in  that,  for  circular  polarization,  in  general  the  peak  spectral  amplitudes 
are  lower  and  the  spectral  widths  are  greater  than  for  vertical  polariza- 
tion. 

Some  of  the  data  runs  exhibit  a second  spectrum  superimposed  on  the 
initial  no.iotonic  continuous  spectrum  of  the  rain.  As  discussed  in  the 
text  (Section  III-B-5)  this  phenomenon  is  caused  by  movement  of  the 
storm  radially  with  respect  to  the  radars.  The  center  frequencies  of 
the  superimposed  spectra  can  be  related  by  the  Doppler  equation  to  the 
wind  speed  recorded  in  the  field  logs  to  within  10Z,  a close  relation- 
ship. Figures  B-5  and  B**29  are  examples  of  this  type  of  spectral  response. 
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Figure  Bl.  Amplitude  of  the  frequency  spectrun  of  a logarithctic 
receiver  for  rain  backscatter;  5 ma/hr  rain  rate, 
9.375  GHz  frequency,  and  W pr larization. 
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Figure  B2.  Anplitude  of  the  frequency  spectrua  of  a logarithmic 
receiver  for  rain  backscatter;  5 m/hr  rain  rate, 

35  GHz  frequency,  and  W polarization. 
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Figure  B7.  Amplitude  of  the  frequency  spectrue  of  e logarithmic 
receiver  for  rain  backscatter;  23  mt/br  rain  rate, 
70  GHz  frequency,  and  W polarization. 
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Figure  B8.  Amplitude  of  the  frequency  apectrua  of  a logarithmic 
receiver  for  rain  backscatter;  23  mU hr  rain  rate, 

95  GHz  frequency,  and  W polarization. 
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Amplitude  of  Che  frequency  spectrum  of  a logarithmic 
receiver  for  rain  backscatCer;  38  tm/hr  rain  rate, 
9.375  GHz  frequency,  and  W polarization. 
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Figure  BIO.  Amplitude  of  the  frequency  spectrum  of  a logarithmic 
receiver  for  rain  backscatter;  38  mt/hr  rain  rate, 

35  GHz  frequency,  and  W polarization. 


Figure  Bll.  Amplitude  of  the  frequency  spectrua  of  a logarithaic 
receiver  for  rain  backscatter;  38  mm/hr  rain  rate, 

70  GHz  frequency,  and  W polarization. 
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Figure  B15.  Anplitude  of  the  frequency  spectreo  of  a logarfth 
receiver  for  rain  backscatter;  70  mm/hr  rain  rate 
70  GHz  frequency,  and  W polarization. 
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Figure  826.  Anplitude  of  the  frequency  spectrua  of  a logarithmic 
receiver  for  rain  backscatter;  70  rm/hr  rain  rate, 

95  GHz  frequency,  and  W polarization. 
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Figure  B17.  Amplitude  of  the  frequency  spectrus  of  a logarichaic 
receiver  for  rain  backscatter;  100  mm/hr  rain  rate, 
9.375  GHz  frequency,  and  W polarization. 
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Figure  B18.  Aaplicude  of  the  frequency  spectrua  of  a logarichaic 
receiver  for  rain  backscatter;  100  tm/hr  rain  rate, 
35  GHz  frequency,  and  TV  polarization. 
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Figure  820,  Amplitude  of  the  frequency  spectrum  of  a logarithmic 
receiver  for  rain  backscatter;  100  nm/hr  rain  rate, 
95  GHz  frequency,  and  W polarization. 
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Figure  B22.  Amplitude  of  Che  frequency  spectrua  of  a logarithmic 
receiver  for  rain  backscatter;  3.1  m/hr  rain  rate, 

35  GHz  frequency,  and  RC  polarization. 


n iu 


m 


JM  MS 


9*7 


739  n* 


li 


Figure  B26.  Amplitude  of  the  frequency  spectrum  of  a logarithmic 
receiver  for  rain  backscatter;  32  na/hr  rain  rate, 

35  GHz  frequency,  and  RC  polarization. 
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Figure  B27,  Aaplitude  of  the  frequency  spectrum  of  a logaritimic 
receiver  for  rain  backs catter;  32  rm/hr  rain  rate, 

70  Cflz  frequency,  and  RC  polarization. 
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Figure  B28.  Aaplitude  of  the  frequency  spectrins  of  a logarithmic 
receiver  for  rain  backscatter;  32  am/hr  rain  rate, 

95  GHz  frequency,  and  RC  polarization. 
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Figure  B31.  Amplitude  of  the  frequency  spectrua  of  a logarithmic 
receiver  for  rain  backscatter;  91  aa/hr  rain  race, 

70  GHz  frequency,  and  RC  polarization. 


